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In this paper we present an analysis of the structure of Bell inequalities, mainly for the case of N 
qubits with two observables each. We show that these inequalities are related to Hadamard matrices 
and define Bell polynomials (in one variable) as an additional tool. With these aids we raise several 
conditions the coefficients of Bell inequalities must satisfy, and recursively generate the whole set 
of inequalities starting from N = 1. Moreover, we prove some characteristic features of this set, 
such as that most of the inequalities contain all expectation values under consideration. Finally, we 
show how the presented results can be used to construct Bell inequalities with certain properties. 
An outlook on further research topics concludes the paper. 



I. INTRODUCTION 

In 1964, J. S. Bell published his now-famous paper pj, 
which revolutionized the foundations of quantum me- 
chanics and our view of nature. After decades of dis- 
cussion Bell demonstrated that it was possible to decide 
experimentally whether the EPR argument Q was cor- 
rect. The central part of this paper was an inequality 
assigned to a specific experimental setup. The expecta- 
tion values of measurements had to satisfy this inequality 
so that a local realistic model could be applied. Even for 
more complex quantum systems there always exists a set 
of linear constraints that serve this purpose. Such con- 
straints are today called Bell inequalities. 

Let us consider two observers, Alice and Bob, who 
perform measurements on qubits (the simplest quantum 
systems, having two possible outcomes). We assume that 
each of them has the choice of measuring one out of three 
observables. According to D. Bohm's variant Q of the 
EPR experiment [2|, these measurements can e.g. be spin 
measurements on particles of spin i , where Alice and Bob 
have the option of measuring at angles 
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respectively (see Fig. QJ. If we assign values to the pos- 
sible outcomes, say +1 to "spin up" and —1 to "spin 
down," we are able to examine expectation values over 
many repetitions of the experiment. Let E\{i) and Ei(j) 
denote the expectation values of measurements at angles 
0i and ijj, respectively, and let E(i, j) denote the expecta- 
tion value of the product of the assigned values. If Alice 
and Bob use an entangled pair of particles from a source 
in singlet state that shows no preference for any specific 
direction, we get 

E l {i)=E i {j) = Q for0<i,j<2. 



Moreover, quantum mechanical calculation |4| yields 



E(i,j) 



if i = j 
otherwise, 



in perfect accordance with empirical data. Thus, if the 
choices made by Alice and Bob are independent of one 
another and each angle is chosen with equal frequency, 
we get the following properties: 

1. If i = J, then the same values are measured. 

2. The overall expectation value of products is zero. 

There is no way of explaining this behavior in terms of a 
local realistic model, where the outcomes of Alice's mea- 
surements are independent of Bob's measurements and 
vice versa. 

Remark. Property [21 is still satisfied if one or both of the 
apparatures are tilted by some arbitrary degree, because 
for any cj> there is 

cos^ + cj)j + cos(7r + (/)) + cos (^" + ^) = 0- 

Thus, essentially the same argumentation holds if 

rji =6i for < i < 2 

is used instead of @, which is often the case in the lit- 
erature (e.g. in 

In a general setup there are N observers having the 
choice of measuring one out of g observables, where the 
outcome of each measurement is q- valued. If we enu- 
merate the observables by 0, 1, . . . ,g — 1, their respective 
choices can be described by a vector 

(fei,fc 2 ,...,fcjv) 

with < ki < g — 1 for 1 < i < N. We can interprete 
this vector as the g-adic expansion of an integer 
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k = (k u k 2) ■ ■ ■ , k N ) g = J^fcj 



N-i 



(3) 
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with < k < g N — 1 . (If there are less than g observables 
at some sites, k simply does not take all possible values.) 
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In a local realistic model each observable is a random 
variable Aj(fcj) in its own right and is independent of the 
choices kj for j ^ i at other sites. We will denote the 
corresponding expectation values of the product of these 
variables by 

E(k) :=E(k 1 ,...,k N ) = (A 1 (k 1 )---A N (k N )y (4) 

which must in any local realistic model satisfy a set of 
Bell inequalities. 

As the systems grow, the number and complexity of 
these inequalities increase dramatically. In fact, in [|| it is 
shown (in terms of joint probabilities, instead of expecta- 
tion values) that the question of whether a local realistic 
model can be applied or not is related to a convex correla- 
tion polytope. The experimental results can be explained 
by a classical probability distribution exactly if the cor- 
responding vector of probabilities and joint probabilities 
lies inside that polytope. This problem is NP-hard. (For 
the definition and a survey of NP-hard problems see pj.) 
Historically, related problems were already investigated 
by G. Boole 8] in the 19 th century. Independently, these 
problems are of relevance in probability theory and re- 
lated research is going on to this day (see also [g] for a 
discussion) . 

Remark. In the literature the enumeration of observables 
usually starts with 1 instead of 0. Therefore, in the lit- 
erature the expectation value @ is written as 



E(h + 1,. 



1) 



with < ki < g — 1 for 1 < i < N. We will also pay 
attention to that convention in this paper, referring to it 
as traditional notation. 



II. BELL INEQUALITIES 

We will now and for the rest of this paper study the 
case of N qubits with two observables each, i.e. q = 2 
and g = 2. For that purpose we consider (see also [Toj]') 
the product 



N 



P(hi, ...,h N ) := Y[(Ai(0) + M*(l)) 



(5) 



for arbitrary hi G {—1,1}, which can be expanded to 

P(h 1 ,...,h N ) = hk i ■■■h h N N A 1 (k 1 )---A N (k N ). (6) 
(k u ... : k N )e{a,i} N 

That again defines a random variable, which now also 
depends on the nonrandom variables hi, . . . , /ij\r. For a 
concrete realization of variables Ai(ki), there is only one 
choice for the hi's so that the product (0 does not vanish, 
in which case each factor is ±2. Thus, we have 



(h 1 ,...,h. JV )e{-i,i} N 



,h N ) = ±2 



N 



Since this sum contains only one nonvanishing term, we 
also get 



,/i«)G{-l,l}' 



,h N )P{hi,. 



IN 



±2 



N 



(7) 



for an arbitrary ±l-valued function c(hi, . . . , /iat). The 
expectation value of that sum must therefore lie between 
— 2 N and 2 N . In order to derive constraints for (0J, we 
substitute 10 in that expression and use linearity of ex- 
pectation. With respect to J2J we set 



(In 



2^ K---h 

.,^«)e{-i,i} JV 



c(hi, h N ) 



(8) 



for < k < 2 



1, and finally get 

a k E{k) 



< 2 



(9) 

By choosing all admissible functions c(hi, . . . , /i_/v) in JSJ, 
the corresponding inequalities 10 represent a complete 
set of Bell inequalities for the experimental setup under 
consideration 0, ^| . That means that these inequali- 
ties are satisfied exactly if a local realistic model can be 
applied. 

Example Let us consider the case N = 2. By (JSJ) we 

have 

a = c(l, 1) + c(l, -1) + c(-l, 1) + c(-l, -1) 
ai = c(l, 1) - c(l, -1) + c(-l, 1) - c(-l, -1) 
a 2 = c(l, 1) + c(l, -1) - c(-l, 1) - c(-l, -1) 
a 3 = c(l, 1) - c(l, -1) - c(-l, 1) + c(-l, -1) 
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with an arbitrary ±l-valued function c(h\,h 2 ). For ex- 
ample, if we choose 



c(hi,h 2 



we get 



(/l! + l)(/l 2 + l) 



■1 ifhi = h 2 = l 
-1 otherwise, 



ao — 2 and a% = a 2 — 03 = —2. 
By (0 this leads, after division by 2, to 

\E(0)-E(1)-E(2)-E(3)\<2, 
which in traditional notation reads as 

\E(1,1) - £(1,2) - £(2,1) - £7(2,2)1 < 2. 



(10) 



(11) 



(The transcription from IjlUI) to Ijlll) happens by writing 
each argument in its binary expansion, using exactly N 
digits, and then incrementing each digit by 1.) By using 
all 2 4 = 16 admissible functions c[h\, h 2 ), we can easily 
verify that this inequality is, up to symmetry, the only 
nontrivial case for N — 2. 

Inequality l|ll|) was first derived b y J . F. Clauser, 
M. A. Home, A. Shimony & R. A. Holt 12], and we will 
therefore subsequently refer to it and its symmetric vari- 
ants as CHSH inequalities. From now on we will also use 
the shorthand notation 



(ao, ai, . . . , a 2 N_ 1 ) 



(12) 



for @. This is convenient since we will see that even 
if J5J is multiplied by an arbitrary nonzero constant we 
can still calculate its upper bound by plain use of (|12fl . 
(Without such multiplication this is trivial, since in that 
case we only need to take the length of this vector to 
achieve this bound.) 



III. HADAMARD MATRICES 

The vector (hi, . . . , Hn) € {—1, 1}^ can also be inter- 
preted as a binary expansion with "digits" ±1. Thus, 
if we use the substitutions 1 t— > and — 1 t— > 1 in this 
expansion, we get a nonnegative integer 



3 = (ii 



,3n) 2 , 



where hi — 1 — 2ji for 1 < i < N. This leads us to 
write Cj instead of c(h\, . . . , hrj), which in the previous 
example means 

c = c(l, 1), ci = c(l, -1), c 2 = c(-l, 1), c 3 = c(-l, -1) 
and further 




We will now study the matrices involved in this opera- 
tion in general. For that purpose we consider the binary 
expansions j — . . . ,Jn) 2 and k = (k\, . . . , kjsi) 2 and 
their scalar product over GF(2), defined as 



N 

(j,k) 2 := ^jih 

z=l 



(mod 2). 



(13) 



(In other words, (j, k) 2 = 1 if the number of l's in which 
the binary expansions of j and k coincide is odd, and 
otherwise.) We observe in (jSJ that = — 1 if and only if 
hi = — 1 and ki = 1, which means that the corresponding 
bits in j and k are both 1. Therefore, we can also write 
© as 



a-k 



i N ~i 

E 

3=0 



with an arbitrary ±l-valued vector (co,ci,. 
Thus, if we define the 2 N x 2 N matrix H 2 n = 



(14) 



• •,C 2 N_ 1 ). 

(hjk) by 



for < j, k < 2 



N 



h jk := (-1)^ 
1, we finally get 



f a ° \ 




( co \ 


ai 




Cl 




= H 2 n 




\a 2 N_ 1 J 




\C 2 N_ 1 J 



(15) 



(16) 



The matrices H 2 n constitute a special type of Hada- 
mard matrices. These are matrices with elements ±1, 
where two rows (and columns) differ in exactly half of 
their elements. (This is another way to say that Hada- 
mard matrices are orthogonal ±l-valued matrices.) They 
have the additional property of having maximal determi- 
nants among all - even complex - matrices with elements 
bound by 1 (in absolute values); namely, if H n denotes 
an arbitrary Hadamard matrix of dimension n x n, then 



|det#„ 



n 2 



(This property comes as no surprise, since this determi- 
nant can be interpreted geometrically as the volume of a 
parallelepiped spanned by the row vectors of H n - each 
of length y^n -, which is maximal if these vectors are or- 
thogonal.) Historically, this was also the reason why such 
matrices were originally studied by J. Hadamard [lJj- It 
is not difficult to prove that for any Hadamard matrix H n 
there is n = 1, n = 2 or n = ik for some k > 1 (e.g. see 
|l4j|): but it is open whether a Hadamard matrix really 
exists for each such n. (The first unknown case up to this 
day is n = 428.) For further information on Hadamard 
matrices see |Q . 

The matrices H 2 n can also be derived 0] via the re- 
cursive definition 



Hi = (1) 



H 2n = 



H n 



for n > 1. 



(17) 
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FIG. 2: The Hadamard matrices H 2 n for 1 < AT < 6. A black square symbolizes +1, a white square — 1. 



(a^) and B, 



This construction was first given by J.J. Sylvester in the 
context of tilings [l6j, 26 years before Hadamard stud- 
ied matrices of that kind. Therefore, these matrices are 
called Sylvester-type Hadamard matrices. In particular, 
they are normalized, i.e. the elements of the first row and 
column are all 1, and symmetric. Using the Kronecker 
(or tensor) product g) for matrices A 
defined as 

'a\\B ai 2 B ■ ■ ■ 
A®B = \ a2lB a22B ' ' ' 



we can write l|17|l for n > 2 also as Bin = H 2 ® B n . 
Therefore we get 



Ho 



Ho 



N 



(The operator (g> is associative, thus it makes no differ- 
ence whether this expression is evaluated from the left or 
from the right.) Explicitly written, that means 




For a visualization of the matrices H 2 n for 1 < N < 6 
see Fig. 

Remark. For matrices of type l|17|) . there also exist nu- 
merous other descriptions; e.g. they can as well be defined 
as the table of characters for elementary Abelian groups 
of order 2 . Moreover, Hadamard matrices have a broad 
spectrum of applications: Hadamard designs [l7| can be 
used for statistical experiments an d also for design- 
ing certain types of bridge tournaments [l9| ; Hadamard 
codes [2(j were used in 1969 by space probe Mariner to 
transmit pictures from Mars to Earth j^J; and Hada- 
mard transforms |22j play an essential role in quantum 
computing, to name just a few. As a consequence natu- 
ral connections arise between the topics mentioned above 
and the results presented in this paper, and vice versa, 
but we will not elaborate on them here. 



IV. GENERAL PROPERTIES OF BELL 
INEQUALITIES 

We will now use certain properties of matrices l|17() to 
derive some general properties of Bell inequalities. First, 
we note that all coefficients ak in @ are restricted to the 
values 0, ±2, ±4, . . . , ±2 , which means, in particular, 
that they are even. Therefore, we can always divide © 
by 2 and still get an inequality with integral coefficients 
and an integral upper bound. Thus, any Bell inequality 
(|9"|) can be written in the form 



]T b k E(k) 



< 2 n 



(18) 



k=0 
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with integral b^s and < n < N — 1. If the b^s are 
relatively prime and 



2"-l 

E 



> 0, 



(19) 



we say that (|18|l is in standard form. (The sum (|19f) can 
never vanish, as we will see below; thus a division of the 
frfc's by their greatest common divisor and an occasional 
multiplication by —1 do the job.) 

We will now derive some general properties of (|18|l . 
regardless of whether it is in standard form or not: 

Lemma 1 The coefficients in 1(18(1 have the following 
properties: 

(i) For each k, there is \bk\ < 2™. 

(ii) // equality holds in (Q) for some k, then bj — for 
all j ^ k. 

(iii) There is always 



2" -1 

E< 

fc=0 



= T 



Proof. We have already noted that the coefficients a k 
in © are restricted to the values 0, ±2, ±4, . . . , ±2^. 
Hereby, the extreme values ±2^ are taken exactly if the 
vector (co, . . . ,c 2 n_ 1 ) in ((16(1 corresponds either to the 
fc-th row vector of H 2 n , or to this vector multiplied by 
— 1. In both cases, we also have a,j = for j ^ k, since 
the row vectors of H 2 n are orthogonal. Division of 
by 2 N ~ n leads to JTSJ, for which now properties (Q) and 
f|TT|) must hold. 

To show (jmj, we denote the column sums of H 2 n by 
£fc for < k < 2 N — 1. By definition of H 2 n we get 

& = 2 N 

£ fc = for k > 0, 



and thus 



E = E = 2 W C 



C . 



(20) 



fe=0 



fe=0 



Taking absolute values in 1(2U|) and dividing by 2 N " 
leads to (fTTTfl . which finishes the proof. □ 



We can easily verify property ljm)l for the CHSH in- 
equality jH) . This property now also justifies the short- 
hand notation 



(bo,h,- 



J 2 N_l 



for (|18l) . since by its use we can always determine the 
upper bound in the corresponding Bell inequality. 

The recursive principle in definition (|17|l enables us to 
take two arbitrary Bell inequalities for N qubits and, by 



using them, construct a new Bell inequality for N + 1 
qubits. For that purpose we define the operator txi as 



(a , . . . , a 2 «-i) X (bo, • ■ ■ , b 2 N_ 1 ) : 
(ao+bo, . . . , a 2 N _ 1 +b 2 N _ 1: ao—bo, 



, a 2 N. 



-6 2 w_i). 



Thus in the first half of the resulting vector the coeffi- 
cients are pairwisely added, and in its second half they 
are pairwisely subtracted. By that we get: 

Theorem 2 Let (ao, . . . , a 2 «_ 1 ) and (bo, ■ ■ ■ , b 2 N_ 1 ) be 
two Bell inequalities for N qubits that satisfy 



2"-l 

E 

k=0 



a k 



2 N -1 

E 1 

fe=0 



Then applying the operator above 
(c , . . . , c 2 jv+i_ 1 ) := (a , . . . , a 2 jv_i) X (6 , 



(21) 



. &2™-l) 



yields a Bell inequality for N + 1 qubits. By substitut- 
ing all Bell inequalities for N qubits in (ao, . . . , a 2 «_ 1 ) 
and (bo, ... , 6 2 w_i) } we get Z?e!/ inequalities for N + 1 

Proof. If we substitute iV i— ► AT + 1 in (|16|l and consider 
the form of the matrix H 2 n+i by setting n = 2 N in 117|l . 
we get exactly the transition described by the ixi operator. 
Therefore the result follows. □ 

By 116(1 we can also show that the rate of Bell inequal- 
ities that do not contain a certain expectation value de- 
creases as N increases. The following proposition quan- 
tifies this fact: 

Proposition The probability that an arbitrarily chosen 
Bell inequality for N qubits does not contain a certain 
expectation value is asymptotically 

1 



as N — > oo. 

Proof. Each coefficient in (0 is by ((16(1 the scalar product 
of a certain row vector of H 2 n and a ±l-valued vector 
(co, . • • ,c 2 w_ 1 ). This product vanishes exactly if these 
two vectors differ in half of their elements. Since the 
number of vectors (cq, . . . , c 2 «_ 1 ) with this property is 



2 N 

2 N-1 



(22) 



division by 2 2 (the number of all inequalities) and Stir- 
ling's formula yield the result. □ 



V. BELL POLYNOMIALS 

It turns out to be convenient to assign to ((18(1 the Bell 
polynomial 



2"-l 



B(z) := ]T b k z k 



(23) 



k=0 
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(These polynomials should not be confused with the mul- 
tivariate polynomials introduced by Eric Temple Bell in 
1934, or the polynomials defined in [llj.) If we want to 
emphasize the number of qubits for which is used, 
we will write B( N '(z) instead of B{z). 

Property (fTTTf) in Lemma ^ tells us that the upper 
bound in the corresponding Bell inequality (|18f) is given 
by |B(1)|; by property Q this is also an upper bound 
for the height of that polynomial (which is defined as its 
maximal coefficient, in absolute values). Multiplication 
of B(z) by an arbitrary nonzero constant does not af- 
fect the corresponding Bell inequality. Therefore we call 
polynomials that can be obtained from one another by 
such a multiplication equivalent. This definition indeed 
leads to an equivalence relation on the set of Bell poly- 
nomials. - Usually we will choose representatives where 
the corresponding Bell inequality is in standard form, 
which means that the coefficients are relatively prime 
and B(l) > 0. (If we do not mind yielding rational coef- 
ficients, we can also consider polynomials with B(l) = 1 
for that purpose.) 

Having taken these preparations, we are now ready to 
adapt Lemma ^and Theorem [21 to Bell polynomials. 

Lemma 3 The coefficients in l|23|) have the following 
properties: 

(i) For each k, there is \bk\ < |£?(1)|. 

(ii) If equality holds in (0) for some k, then bj = for 
all j 7^ k. 

Proof. The lemma follows immediately from the defini- 
tion of Bell polynomials and Lemma ^ □ 

For the adaptation of Theorem [5] we observe that 

(a , a fc _i, 0, 0) <-> A(z) 

(0, 0, a , ojfe-i) <-> z k A(z). 

This leads us to define 

A(z) tx B(z) := (1 + z 2 ™) A(z) + (l - z 2 ")£?(z), 

where A(z) and B(z) are Bell polynomials for N qubits. 
By that we get: 

Theorem 4 Let A{z) and B{z) be two Bell polynomials 
for N qubits that satisfy 

|A(1)| = |B(1)|. (24) 

Then applying the operator above 

C(z) := A(z)txiB{z) 

yields a Bell polynomial for N + 1 qubits. By substituting 
all Bell polynomials for N qubits in A(z) and B(z), we 
get all Bell polynomials for N + 1 qubits. 



Proof. The theorem follows immediately from the defini- 
tion of Bell polynomials and Theorem [21 O 

Example Let us consider the two CHSH inequalities 

(1,1,1,-1) and (1,-1,-1,-1). (25) 

They obviously satisfy l|21|l . hence we can apply Theo- 
rem^ and get 

(1,1,1,-1)^(1,-1,-1,-1) = (2,0,0,-2,0,2,2,0). 

After division by 2, in traditional notation this reads as 

\E(1, 1, 1) - E(l, 2, 2) + E(2, 1, 2) + E(2, 2, 1)| < 2, 

which is an MABK inequality [MHHI. Alternatively, 
we can use the corresponding Bell polynomials and The- 
orem 0] to achieve the same result. In that case we have 

(1, 1, 1, -1) <-> A(z) = 1 + z + z 2 - z 3 
(1, -1, -1, -1) «-> B(z) = 1 - z - z 2 - z 3 

and 

A(z) xi B{z) = (1 + z 4 )(l + z + z 2 - z 3 )+ 
(1 - z 4 )(l - z - z 2 - z 3 ) 
= 2-2z 3 + 2z 5 + 2z 6 , 

which indeed corresponds to the inequality derived 
above. (Note also that whenever l|21|l is satisfied, the 
same is automatically true for the corresponding Bell 
polynomials and l|24|).) 

If we replace the second CHSH inequality in 1251) by 
the trivial inequality (1,0,0,0), we must first multiply 
this inequality by 2 in order to satisfy condition 121fl . 
Application of Theorem [21 now yields 

(1,1, 1,-1) 1*1(2, 0,0,0) = (3, 1,1, -1,-1, 1,1,-1), 

which corresponds to 

3 + z + z 2 - z 3 - z 4 + z 5 + z 6 - z 7 . 

(This time we leave it to the reader to obtain that poly- 
nomial directly by use of Theorem^) In traditional no- 
tation, this reads as 

|3£(1, 1, 1) + E(l, 1, 2) + E(l, 2, 1) - £7(1, 2, 2) 

- £7(2, 1, 1) + £7(2, 1,2) + £7(2, 2, 1) - £7(2, 2,2)] < 4, 

which is already in standard form. 

We will now study the structure of Bell polynomials 
in detail. For N = 1, the whole set of polynomials is 
given by ±1 and ±z. Theorem 01 tells us that the Bell 
polynomials for N — 2 have the form 

±1 ±1 

(l + z 2 / + (l-z 2 / 

\ \ 
±z ±z 
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(presented in a hopefully self-explanatory notation) . Pro- 
ceeding in that way, we see that B( N > (z) consists of the 
2 JV_1 summands 



Theorem 5 Let sj^ (z) be as defined in (|28|l . Then the 
complete set of Bell polynomials for N qubits is given by 



(26) 



with each of them multiplied by a factor ±1 or ±z. 
Since there are exactly four choices for this multiplication 
within each summand, the number of Bell polynomials is 
indeed 4 2 = 2 2 . 

It turns out to be convenient to choose a fixed enu- 
meration for the summands (|26|) . Therefore, we define 

8?\z) := (1 + (-l) fe "- 2 z 2 "") • • ■ (1 + (-l) k "z 2 ) (27) 
with 



(k 



N-2, ■ 



,ki,k, 



0)2 



N-2 
i=0 



whereby the empty product is set 1 as usual. Thus, we 
have 



(31) 



fc=0 



where u = (u 2 n-i_i, . . . , 1*0)2 and v = . . . , ^0)2 

are arbitrary numbers with binary expansions of length 
2 N ~ 1 (occasionally written with leading zeros). 

Proof. We have already achieved the structure of these 
polynomials in the recursive process that led to the def- 
inition of s^\z). The different choices of factors ±1 
and ±z at the k-th summand are now described by the 
Boolean variables u k and v k , so we are done. □ 

Remark. A more direct approach would be to define the 
polynomials 

t[ N \z) := h ok + h lk z + h 2k z 2 + ■■■ + hx-i^z"- 1 (32) 



s (1) 


[z) 


= 1 








(2) 
s 


[z) 


= H 


-z 2 








[z) 


= 1 - 


-z 2 








[z) 


= (1 


+ ^ 4 )(1 








[z) 


= (1 


+ ^ 4 )(1 






s (3) 


[z) 


= (1 


-z 4 )(l 








[z) 


= (1 


-z 4 )(l 






4 4) 


[z) 


= (1 


+ z 8 )(l 


+^ 4 )(1 + 


J) 


4 4) 


[z) 


= (1 


+ z 8 )(l 


+^ 4 )(1- 


z 2 ) 



We can also use a recursive definition for s^\z), namely 
4\z) :=1 

{N) _Ul + z* N - 1 )s[ N - 1 \z) iffc<2^ 2 
[{l-z A )4V-2( Z ) if A; > 2 iv 

with < k < 2 7V ~ 1 . Expanding (z) yields an expres- 
sion of the form 



l±z 2 ±z 4 ±---±z 2 - 2 , 



(29) 



thus si N \z) is an even polynomial with coefficients ±1 
and degree 2 N — 2. Multiplication by z yields 



z ± z 3 ± z 5 ± • • • ± z 2 



by using 1|15[) . and consider 

Bi N \z) -JeVi) 



AN) 



fc=0 



for an arbitrary u; = (it^jv^, 100)2- However, the 
form (|31|l reveals more of the structure of these poly- 
nomials, which turns out to be useful. The connection 
between (|2Tjl and lp!2^l is given by 



,(A0 



(33) 



for < k < 2 N -\ 



Example For N — 1 and N — 2 the complete set of Bell 
polynomials is given by 



= 1 

B«(z) = z 



< = -1 
B® = -z 



and 



B$(z 
B$(z 
B$(z 

B{ 2 o\z 
B[l\z 
B$(z 



(30) B[ 3 '(z 



2 

1 + z 
1 + z + z 2 
2z 

-2z 2 

1 - z - z 2 
-1 + z - z 2 
-2z 3 



B 



z 2 + z 3 



B${z 



B. 



22 \ z 



Bil'iz 

Bil\z 
Bgiz 



2z 2 

-1 + z + z 2 + z 
1 - z + z 2 + z 
2z 3 

-2 

— 1 — z + z 2 — z 
-1 - z 
-2z. 



z 2 + z 3 



which now is an odd polynomial with coefficients ±1 and 
degree 2^—1. 



Note again that all nontrivial cases above correspond to 
CHSH inequalities. 
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Lemma 6 Let \z) be as defined in (|31[1 . The 

€'(1) =(-l)«o 2 ^-l 
^(-1) = (-l)"o+«o 2 W-l 



fc=0 



Furthermore we have 

B {N) (-z\- B (N) (z) 

where u and u® v are defined as 

u := (1 - u 2 «-i-i 5 ■ • 1 - "0)2 
w © v := (m 2 «-i-i © W2«-i-i, ■ ■ ■ , uo © ^0)2 



(34) 
(35) 

(36) 



(37) 
(38) 



(39) 
(40) 



and the operator © on the right-hand side of 14U|) means 
addition over GF(2). (Thus, u(Bv defines a bitwise "ex- 
clusive or" operation on u and v.) 



Proof. For all k > there is 



(1) 



(-i) = o, 



since in that case si (z) always contains a factor 1 - 
for some j. Therefore, by (|31|l and 



(41) 



4 N \l) = 2 N - 1 

we immediately get Q34JI and (|35|l . Similarly, for any A'' 
and k there is 

4^(0) = 1. 

Thus, by setting z — in (|3 1|) we get the following: The 
fc-th summand contributes only to the sum if Vk = 0, 
namely +1 if = and —1 if = 1. This is exactly 
what (|36[1 tells us formally. 

The transformation B^\z) > — B^\z) can be 
achieved by substituting Uk >— » 1 — Mfe for < k < 2 Ar ~ 1 , 
since the latter changes the sign of every summand in 
(|31ll : therefore 13711 holds. Finally, the transformation 
z 1— > — z changes the sign of the fc-th summand exactly 
if Wfc = 1, since s^(z) is even for any N and k. This 
shows that (I38|) holds, and thus completes the proof. □ 

Corollary The polynomials (|31|l satisfy 

|5W(1)| = |SW(-1)|=2 JV - 1 . (42) 

Furthermore, B^ (z) is even exactly if v — 0, and odd 
exactly if v — 2 2 — 1 . 

Proof. Property (|42|l is a trivial consequence of i|34[l and 
is even, then by (|38|) we get 

B( N) (z)-B [N) (z) 



Hence u = u © v, which means that 

u = (0,0,...,0)a = 0. 

Conversely, if v = 0, then B^\z) = B { u ^\-z) by J3HJ|, 
and Bu^\z) is thus even. 

On the other hand, if B^\z) is odd, then by (|37|l and 
we get 



B (N) (z) - B ( "> (z) 
Hence u = u © w, which means that 
u = (l,l,...,l) 2 =2 2 



1. 



Conversely, if t> = 2 2N 1 — 1, then B^\—z) = —B^\z) 
by (|3"T|) and ||3SJ, and B^,\z) is thus odd. So we are 



done. 



□ 



Because of (|42() . it is sometimes useful to consider the 
normalized Bell polynomials 

^W^^CW, (43) 
which have the property that 

|5W(1)| = |5W(-1)| = 1. 
For 1 < N < 3, they are depicted in Fig. 



VI. ANALYZING THE STRUCTURE OF BELL 
INEQUALITIES 

The numbers u and v completely determine the poly- 
nomial (|31|l for any fixed N. Since the bits in the binary 
expansion of u coincide with the signs of the summands 
in i|31|) . we call u the sign number of Bu^\z). On the 
other hand, the bits in the binary expansion of v describe 
whether a summand in l|31|l is an even or an odd poly- 
nomial (of form (|29|l or (|30|l . respectively). Therefore we 
call v the parity number of B^\z). 

Since B^\z) and B^\z) are equivalent by (|3TJl . it 
is sufficient to consider even sign numbers. (This is par- 
ticularly the case if the corresponding Bell inequality is 
in standard form, since from B$P(1) > and it fol- 
lows that uq = 0.) Furthermore, we can extend the no- 
tion of equivalence to enclose symmetry transformations 
as well. This means that we do not distinguish between 
inequalities that can be transformed into one another by 
permutations of sites, observables or measurement values 
(since those can be considered as not "essentially differ- 
ent"). This clearly reduces the number of equivalence 
classes for any fixed N, though it does not substantially 
affect the growth rate of this number as N — > 00. 

We will now study properties of Bell polynomials in 
terms of u and v. But before we do this we make some 
simple observations: 



-1 -0.5 0.5 1 



-1 -0.5 0.5 1 



-1 -0.5 0.5 1 



FIG. 3: The normalized Bell polynomials B ( U J (z) for 1 < N < 3. 



Lemma 7 Let s k (z) be as defined in l)28[l. Then: 

(i) All coefficients in s N \z) are +1. 

(ii) For k > half of the coefficients in s^\z) are +1, 
the other half being — 1. 

(iii) For each N , there is 

E ^\z) = 2^-\ 

fc=0 

Proof. By using , the content of this lemma is an im- 
mediate consequence of the properties of the Hadamard 
matrices investigated in Section ITY1 We can also prove 
this lemma directly: Property (0) is obvious since s N ^ (z) 
consists only of factors 1 + z 2,1 , which all have positive 
signs. By |{4*T)) we have already seen that s < k N) (l) = for 
k > 0. If, on the other hand, we set z = 1 in (|29|l . we get 
an expression of the form 

1±1±1±---±1. 

This sum can only be zero if the number of +l's equals 
the number of — l's, thus |p} holds. Finally, (fTTTfl can be 
proved by induction (which would have also been possible 
in case of lpi]l). Obviously lfm)l holds for N = 1, so let us 
assume by induction hypothesis that it is true for N — 1. 
By setting 

Sn:= £ 4*>(*), 

fc=0 

this means that Sn-i — 2 JV ~ 2 . Now by we get 
= (l + z 2 )5 JV _ 1 + (l-z 2 )S N .x 

no r)iV-l 
= Zjpf—l = Z 

This concludes the proof. □ 

Hint. Property (JmJ can be generalized for arbitrary sums 
of terms 

(l±z fel )(l±z fc2 )---(l±z fe ") 



over all possible sign combinations. By symmetry we see 
that this sum is always 2", independently of the concrete 
values of the fcj 's. We can therefore immediately simplify 
an expression like 

(1 + z 8 )(l + z 2 ) + (1 + z 8 )(l - z 2 ) + (1 - z 8 )(l - z 2 ) 

to 

4-(l-z 8 )(l + z 2 ) = 3-z 2 + z 8 + 2 10 , 

since the sum over all terms (1 ± z 8 )(l ± z 2 ) must be 4. 
The sketched method can frequently be applied in the 
process of determining Bell polynomials (see the example 
below) . 

Theorem 8 For any Bell polynomial 

2 N -l 

B uv (z) = (44) 

fe=0 

as defined in (|31|l . we have: 
(i) If u — 0, then 

b = - h 

h = -&3 
h = -h 



b 2 N_ 2 — —62^-1- 

(ii) If v = 0, then b 2 j+i = for all j > 0. 

(iii) If v is even, then 

b\ + 63 H 1- b 2 N_ 1 = 0. 

(iv) If v is odd, then 

bo + h-\ h b 2 n_ 2 = 0. 

(v) If the binary expansion of v contains an even (odd) 
number of 1 's, all bk 's are even (odd). 
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Proof. Let 

2 Jy - 1 -l 

p(z):= (-1)^(1-^)4^^) (45) 

2 N - 1 -1 

q(z):= (-lT k v k si N) (z), (46) 

then 

B uv (z) = + zq(z). 

In other words, p(z) denotes the even part of B uv (z), and 
zq(z) denotes its odd part (because p(z) and q(z) are 
both even). If u = 0, then by property ljm|) in Lemma [7| 
we have 

p(z) = - q(z), 

therefore Q holds. On the other hand, if v = then 
q(z) = 0, and now (JijJ follows (since in that case the 
polynomial B uv (z) is even). 

To show (JmJ and O we use (|34|) and (|35|) . by which 
we get 

bo + h + b 2 +- ■ ■ + b 2N _ 2 + 6 2 «-i = {-l) u °2 N ~ l 
b -W+ b 2 - ■ ■ + b 2N _ 2 - b 2N _, = (_i)«o+«o 2 JV-^ 

Subtraction and addition of these equalities now prove 
(|iu)l and l(iv|) , respectively (since vq = if v is even, and 
vo = 1 otherwise). Finally we get @ by observing that 
an even (odd) number of nonvanishing summands in l|45(l 
and (|46[) also leads to even (odd) coefficients in these 
polynomials. □ 

Remark. The theorem above can be enhanced in the 
following ways: 

(a) Properties iJuJ-lO are also valid for equivalent poly- 
nomials of B uv (z), whereas properties and (0) 
need some slight adaptation if equivalent polynomials 
of B uv (z) are considered. 

(b) Similar results to (0) and (JuJ can be derived for other 
specific values of u and v, such as 

u = 2 - 1 or v = 2 2 - 1. 

(c) If B^J (z) denotes (by a short-term overload of nota- 
tion) the fc-th derivative of B uv (z), then 

h = ~Bg>(0). 

We can therefore formulate properties of the b^s also 
in terms of derivatives of B MV (z); e.g. property 10) 
then reads as 

B' Qv (0) = 2 N - 1 -B Ov (0) 
B^ +1) (0) = -(2j + l)4f(0) for j > 1. 



(d) The frfc's can also be written in terms of u and v. For 
instance 

b = (u,v) - (u,v), (47) 

where (u, v) denotes the scalar product of the binary 
expansions of u and v (which counts the positions of 
l's at which the binary expansions of u and v coin- 
cide). 

(e) By l|47l) we can again count the Bell inequalities that 
do not contain E(0), as was done for arbitrary ex- 
pectation values in the proof of the proposition at 
the end of Section llVl Now bo — holds if and only 
if exactly half of the l's in the binary expansion of 
u coincide with l's in the binary expansion of v] to- 
gether with (|22[1 , this proves the formula 

Example Let N — 4 and 

u = (0,0,0,0,1,1,1,0) 2 
v = (0,0,0,0,0,0,0,0) 2 . 

Then 

B uv {z) = (l + z 8 )(l + z 4 )(l + z 2 ) 

-(l + z 8 )(4-(l + z 4 )(l + z 2 )) 
+ 4(l-z 8 ) 

= 2(1 + z 2 + z 4 + z 6 - 3z 8 + z 1Q + z 12 + z 14 ), 

using the hint provided after the proof of Lemma We 
can now easily verify properties (JnJ and iQ) of the theo- 
rem above. (Actually, property JmJ also holds, but this is 
trivial if (JnJ can be applied.) If we consider alternatively 

u = (0,0,0,0,0,0,0,0)2 
V = (0,0,0,0,1,1,1,0)2, 

we get 

B uv (z) = (l + z s )(l + z 4 )(l + z 2 ) 

+ z(l + z 8 )(4-(l + z 4 )(l + z 2 )) 
+ 4(l-z 8 ) 
= 5 + 3z + z 2 - z 3 + z 4 - z 5 + z 6 - z 7 
- 3z 8 + 3z 9 + z w - z 11 + z 12 - z 13 + z 14 - z 15 . 

This time properties JJJ, Jm)l and (Q| of the theorem 
above can be verified. 

From now on we will call a Bell inequality full-term if it 
contains all expectation values under consideration. (We 
can further call a Bell inequality t-term if it contains pre- 
cisely t expectation values; then the full-term inequalities 
for N qubits are 2 JV -term, and the trivial inequalities are 
exactly the 1-term inequalities.) Using that diction we 
get: 
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Theorem 9 Any complete set of inequivalent Bell in- 
equalities for N qubits has the following characteristics: 

(A) Exactly 2 N inequalities are trivial. 

(B) At least half of the inequalities are full-term. 

Proof. The number of expectation values that can pos- 
sibly appear in a Bell inequality for N qubits is 2 . Up 
to equivalence, exactly one trivial Bell inequality corre- 
sponds to any of these values; therefore (0) follows. 

By property (J]J in Theorem |H1 we know that an odd 
number of l's in the binary expansion of v leads to odd 
coefficients in Q44|l. In particular, these coefficients are 
therefore all different from zero. Since half of all v's have 
the property stated above, at least half of the correspond- 
ing Bell inequalities are full-term. This fact still holds 
if equivalent polynomials are eliminated from that set, 
which proves JbJ. □ 

Actually, for N — 2 and N = 3 exactly half of the Bell 
inequalities are full-term. Thus, in these cases there is 
a one-to-one correspondence to v's with an odd number 
of l's in their binary expansion. It follows that each of 
these inequalities, in standard form, has the upper bound 

We give a concluding example in order to illustrate 
how the results above can be used to construct Bell in- 
equalities with certain properties: 

Example Find all Bell inequalities for N qubits in stan- 
dard form, where the coefficient of E(k, k, . . . , k) is max- 
imal. 

Solution. We have already listed the Bell polynomials for 
N = 1 and N = 2, so let us assume that N > 3. We will 
first consider the case k = 0. By (|T%|l and property 10) 
in Lemma n we know that the upper bound in any such 
inequality, in standard form, is less than or equal 2 N ~ 1 . 
Property (uj in the same lemma tells us that the coeffi- 
cient bo = 2 N ~ 4 can hereby never be admitted, since in 
that case the inequality can always be reduced to 

|£(0,0,...,0)| < 1. 

Thus we are aimed to consider bo = 2 N ~ 1 — 1. 

By property (0) in Theorem [S] all coefficients of the 
corresponding inequality are odd; otherwise, this inequal- 
ity could again be reduced and property (Q) in Lemma ^ 
would be violated. Consequently, all inequalities of the 
desired type are full-term. The former property tells us 
further that the binary expansion of v contains an odd 
number of l's. Together with (|47|l . which reads as 

(u,v)-(u,v)=2 N - 1 -h 

this leads to 

(u, v) = 2 N - X - 1 and (u, v) = 0. 

Recalling the definition of standard form, the correspond- 
ing Bell polynomial B(z) must also satisfy B(l) > 0; 



thus, by l|34|) we further have uo — 0. All in all, that 
means the solutions for k = correspond to 

u= (0,0,..., 0,0,0)2 and v = (0,0,..., 0,0,1) 2 , 
u= (0,0,...,0,*,0) 2 and v = (0, 0, . . . , 0, 1, 0) 2 , 
u= (0,0,...,*, 0,0)2 and u = (0,0,..., 1,0, 0)2, 



u= (0, *,..., 0,0, 0)2 and v = (0, 1, . . . , 0,0, 0) 2 , 
u= (*,0,..., 0,0,0) 2 and w = (1,0,..., 0,0, 0)2, 

where * denotes an arbitrary bit. Hence the number 
of corresponding Bell polynomials is 2 N — 1. The poly- 
nomials for k = 1 are finally reached by the symmetry 
transformation 

B(z)^z 2N - 1 b{^), 

since this is exactly what happens if the enumeration of 
observables is reversed (and thus k 1— * 1 — k). 

To see a concrete instance of polynomials, let us con- 
sider the case N = 3. The desired Bell polynomials for 
k = are 

3 + z - z 2 + z 3 - z 4 + z 5 - z 6 + z 7 
3 + z + z 2 - z 3 - z A + z 5 + z 6 - z 7 
3 - z + z 2 + z 3 - z 4 - z 5 + z 6 + z 7 
3 + z- z 2 + z 3 + z 4 - z 5 + z 6 - z 7 
3 - z - z 2 - z 3 + z 4 + z 5 + z 6 + z 7 
3 + z + z 2 - z 3 + z 4 - z 5 - z 6 + z 7 
3 - z + z 2 + z 3 + z 4 + z 5 - z 6 - z 7 , 

and the remaining polynomials for k — 1 are yielded by 
the transformation 

B(z) h+ z 7 B^~y 

(We leave it to the reader to write down the correspond- 
ing Bell inequalities in traditional notation.) 

If N is small, we can solve analogous problems also 
by brute force. But even for not too large N, this is 
hopeless since the number of Bell inequalities not only 
grows exponentially, but swperexponentially. With the 
tools provided in this paper, however, specific results can 
even be obtained for large N. 

VII. OUTLOOK 

So far nothing has been said about quantum violations, 
which can also be investigated in the presented context. 
Hereby, the maximal violation is always obtained for a 
generalized GHZ state 0, |2(| . We also mentioned that 
by taking symmetry transformations into consideration, 
the set of "essentially different" Bell polynomials is re- 
duced; so a closer look at the structure of this condensed 
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set is interesting. And, of course, we can also study more 
general setups (and the corresponding quantum viola- 
tions), where each measurement has more than two pos- 
sible outcomes or each observer has the choice of more 
than two observables. 

Apart from being of theoretical interest for the founda- 
tions of quantum mechanics, Bell inequalities are also of 
practical use in quantum cryptography (which is some- 
times more accurately termed quantum key distribution) . 
A. K. Ekert j2]j presented a variant of the method devel- 
oped by C.H. Bennett & G. Brassard |2l|. according to 
which an eavesdropper may be recognized by checking 
whether a certain Bell inequality holds. This approach 



was recently expanded by V. Scarani & N. Gisin to quan- 
tum communication between N partners pflj . where the 
security of this communication is again linked to viola- 
tion of Bell inequalities. 
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